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Abstract 

We extend the discussion of the growth kinetics of the large-N 
time-dependent Ginzburg-Landau model with an order parameter 
described by a <I> 6 free energy functional, to the conserved case. 
Quenches from a high temperature initial state to a zero tem- 
perature state are studied for different selections of parameters 
entering the effective potential. In all cases we obtain the asymp- 
totic form of the structure factor. As expected in analogy with 
the well studied $ 4 model, we find multiscaling behavior whenever 
stable equilibrium is reached in the ordering region. On the other 
hand the present model also displays a novel feature, namely the 
occurrence of metastable relaxation. 



Introduction 

In a previous paper M, hereafter referred as I, we examined the 
growth kinetics of an N-component field with a non-linear local interac- 
tion of the sixth order, namely a $ 6 term, subject to a sudden quench 
from a high temperature initial state. The equations of the large N spher- 
ical model have been deduced and solved in the case of a non-conserved 
order parameter (NCOP or model A) at zero temperature. 

In the present paper we turn our attention to the kinetics of the 
conserved parameter case (COP or model B). 

In contrast with the non conserved case one can observe multiscal- 
ing behaviour whenever stable equilibrium configurations are reached 
in the presence of a potential with two degenerate minima (see I, fig.l). 
The breakdown of standard scaling || is due to the same mechanism as 
for the $ 4 model ||: two distinct dominant lengths compete in the late 
stage of the dynamical process, namely the dimensional length L(t) = 
(2rt) 1 / 2 , with z = 4, and the wavelength A;" 1 of the peak of the structure 
factor. Since these two lengths differ only by a logarithmic factor one 
cannot speak of a single dominant length. As a consequence, each mode 
k evolves with a different exponent, a behaviour which has been named 
multiscaling. 

In analogy with the NCOP, the <3> 6 model exhibits also metastable 
solutions which are absent in the spherical $ 4 . In the following section, 
in fact, we shall extend a theorem concerning the metastability of the 



model which is the equivalent of the one already proved for the NCOP. 
This suggests that the metastability is not due to the kind of relaxation 
dynamics assumed, but it is a property of the functional form of the free 
energy. 

The presence of metastability is interesting because, as stressed in 
I, it is commonly believed that the metastable solutions disappear when 
taking the limit iV — ► oo. 

The present paper is organized as follows: 

In section 1 we solve the model by explicit integration of the equations 
of motion for the structure factor, in the asymptotic regime. Since these 
equations have been obtained in the first part of this work we refer to I 
for a more detailed presentation of the model and of its large- N limit, as 
well as for a discussion of the initial conditions. Finally in the last part 
we summarize, discuss some result and conclude the work. 

1 Solution of COP the model 

We start by considering the conserved model (p = 2) explicitly and solv- 
ing the equation of motion for the structure factor, eq.(12) of I. In the 

_ -x 4 

present case, introducing the scaling function F(x) = ^—g- of the dimen- 
sionless wave- vector x = kL(t), eq. (10) of I can be formally integrated, 
yielding: 

C(k, t) = Ae~ mx2 L e (t) F(x) (1) 



(2) 



where: 

P(t) = sign(Q(t)) 

is the squared rate between the two lengths C(t) = [2r|Q(t)|] 1//2 and L(t), 
and Q(t) is obtained by self-consistency eliminating S(t) from eq.(9) with 



m 

L(t) 



the help of eq.(10) of I. From this equation we obtain: 



S(t) = AL(t) e - d fd d xx- e e- 



(3(t)x 2 -x A 



(3) 



In order to calculate the integral on the right hand side one has to con- 
sider three different cases, according to the behaviour of (3(t) in the late 
stages (3|: 

a) limt-Kx, /3(£) = 



In this case L(t) prevails over C{t). To first order in j3(t), one obtains 
asymptotically: 

S(t) ~ AL(tf- d f d d xx- e [l - f3(t)x 2 ]e~ xi = AL(t) e - d [I - (3{t)h] (4) 



where: 



and: 



I = f d d x x~ e e 



d„ Jl-6 „-ai 4 



h = / d d x iT v e 



(5) 



(6) 



b) lim i ^ 00 /5(t) = oo 



Now C(t) is asymptotically dominant. Since the term j3{t)x 2 prevails 
over x 4 in the exponential of eq. (H), one can disregard the latter and 
obtain: 

S(t) ~ AL(t) - d f d d x x- 9 e~ mx2 ~ AL(t) 9 - d (3(t)^ (7) 

in the late stages of the dynamics. 
Finally one considers: 

c) lim^oo (3(t) = -oo 

In this case S(t) can be evaluated asymptotically by the steepest 
descent technique: 

S{t)rs,L{t) 9 - d Ae^^t) e -^ (8) 

Because of this threefold possibility in the calculation of S(t) we need 
to know the qualitative behaviour of S(t) in the late time regime. As 
in the non conserved case [1], it is sometimes (e.g. in cases /X3 and fi^) 
not a priori evident whether the relaxation into the metastable $ = 
configuration is to be expected. Here we are able, as for NCOP, to 
produce a criterion which establishes under which conditions metastable 
solutions are allowed. 

A necessary and sufficient condition in order to observe relaxation in 
the $ = final state (i.e. lim^oo S(t) =0), is to have: 

d 2 V(<P) . 



\4>=o 



>0 (9) 



and: 



m < -s-^^f ) (10) 



where t is a generic time instant in the asymptotic region. 

We outline the proof of the criterion beginning with the necessary 
condition eq. (H), i.e: 

{hmS(t)=0}^{^i|^ >0} (11) 

In fact, since: 

d 2 VU) . , N 

-^U-^r (12) 

in the case r < 0, recalling that lim^oo S(t) = 0, eq.(9) of I reads: 

Q(t)~r (13) 

for long times. From eq. (0) one sees that limt_ i . 00 /3(t) = — oo. Inserting 
this result into eq. (§), we obtain: 

S(t) ~ e^Lit) 2 ^ (14) 

which is not consistent with the assumption lim^oo S(t) = 0. We con- 
clude that it must be r > 0. The explicit solution of the model proves 
that consistent solutions effectively exist in this case (see section 2). 

Secondly, as regards sufficiency, both eqs. @ and fllPf) are required 
simultaneously. In fact, from eq.(9) of I we deduce that, if the condi- 
tion (Hg) is fulfilled, then: 

Q(t) > (15) 

6 



We distinguish three asymptotic scenarios: 

eL)lvoat-+oo P(t) = 0. 

In this case, from eq. @), lim^oo S(t) = 0. [q.e.d.] 

b)Um t _, tX3 (3(t) = oo. 

Now, from eq. (|^), we have: 

S(t) _ d-9Q(t) 

S(t) " l 2 >Q(t) l ° J 

In this case, since S(t) is positive defined, from eq. ([15|) and eq. (|T6| ) we 
deduce that S(t) < and therefore lim^oo S(t) = 0. [q.e.d.] 

cjlim^oo /3(t) = -oo. 

In this case, from eq. (||) we get: 

S{t)_ d-6 Q(t) 2TQ\t) 1 \Q(t)\ 

S(t) l 2 >\Q(t)\ [ d-6 U+ 2t Q( f ) j + iJ {U> 

In this case, since lim^oo Q(t) = — oo from eq. ( p~5| ) and eq. ( |TTD we 
deduce that S'(t) < and so lim^oo S(t) = 0. [q.e.d] 

To summarize we observe that also in this case, as for NCOP, the 
possibility of metastable relaxation, established by this criterion, is due 
to a local property of the functional form of the potential around the 
metastable solution, i.e. eq. (H), and to a dynamical property of the 



asymptotic regime, i.e. eq. fllPp. In practice a numerical solution of 
eq. (12) of I reveals that, when eq. (|S]) is fulfilled, it is always possible 
to achieve the condition ( |TUD by decreasing the variance of the initial 
condition 5(0). 

We turn now to the solution of the model by considering different 
cases, according to the parameters \i = (r, g, A) characterizing V($) (see 
I, fig.l). 

Let us begin with the case of simple diffusion: 

Ho = (r = 0, g = 0, A = 0) 

This case is trivial to compute since Q(t) =0, but interesting because 
the existence of a fixed point at /io affects the behaviour of dynamical 
processes characterized by different values of the parameters /x. By ex- 
plicitly calculating the structure factor, from eq. ([!]) we find: 

C{k,t) = AL 9 {t)F{x) (18) 

This result, which is exact at all times (and for all N), shows that scaling 
holds true from beginning to end, as in the non-conserved case. 

fii = (r = 0,g = 0,X>0) 

This choice of the parameters represent the tricritical || case which 



can be solved by considering, firstly, the case a) lim t _ i . 00 (3{t) = 0. From 
eq.(9) of I and eq. (H), to leading order, one obtains: 



Q(t) ~ aL(ty (e - d+2) + c (19) 



where a = T , d _ d ° i) and c are constants. 



From eq. (|i~9"|) we evaluate: 

P(t) ~ L{tf {e - d+l) + cL~ 2 (t) (20) 

Therefore this solution is consistent only for d > d c , with d c — 9 + 1, as 
for NCOP. In this case we find: 

C(k, t) ~ Ae-[ aL W 2( " c " d)+cL W^ 2 L 9 (t) F(x) (21) 

Hence: 

C(jfe, t) ~ A[l - aL(t) 2 ^~ d) x 2 }L\t) F(x) (22) 

for d c < d < d c + 2, while: 

C(jfe, t) ~ A[l - 5L(t)- 2 x 2 ]L e (t) F(f) (23) 

for d > d c + 2 (where c = c for d > d c + 2). 

So, for d > <i c , the trivial fixed point at /x is still attractive and 
the system scales asymptotically as in the /i case, with x dependent 
corrections to scaling. 

We consider now quenches at d < d c . In this case lim^oo (3(t) = oo 
since (3{t) — > — oo is not allowed for fii. From eq.(9) of I and eq. (|7|) we 
deduce: 

Q(t)~[(d-6+l)(aL 4 (t) + c)}^tri (24) 
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with a = , r\ A -? , „ and c is a constant. Therefore: 

(2r) d ~ d c+2 

P(t) ~ L(£) 2 ^Ti (25) 

which is consistent for d < d c . In this case C(t) prevails asymptotically 
over L(t) and we have to look for scaling with respect to the former 
length. Therefore we go back to eq. (1), which can be written as: 

C(k, t) = Ae- k4L4{t) C e (t)F(S:) (26) 

where the scaling function is now defined as: 

n%) = ^ (27) 



:r 



and the dimensionless variable x = kC(t) is expressed in terms of the 
dominant length C(t). Substituting, with the help of eq. (|25l), L(t) with 
C{t) in eq. ( ^6|) we obtain: 



C(k,t) = Ae- aC2(d -^C 9 (t)F(x) 



A 



1 - aC 2(d -^ )i4 ] C\t)F{x) (28) 



where a is a constant. This result shows that scaling holds true even 
in this case, with x dependent corrections, but with a modified scaling 
function. The power growth law of the dominant length, C(t) ~ £i, is 
obeyed with an exponent z — 2{d — d c + 2), which depends on the space 
dimensionality. 

// 2 = (r = 0, g > 0, A > 0) 
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In this case, as for NCOP, we expect the parameter A to be irrelevant 
and the asymptotic form of the structure factor to be the same as for a 
$ 4 theory (with r = 0) (see ||). In fact, since lim^oo S(t) = 0, solving 
eq. (9) of I for long times: 

Q(t)~gS{t) (29) 

we find: 

C(k, t) ~ Ae~ 9Aa ' m :+ x2 L e (t) F{x) (30) 

for d < d c = 6 + 2 

C{k, t) ~ A[l - gAaL{t) d *- d x 2 }L e {t) F{x) (31) 

for d c < d < d c + 2, and: 

C(jfe, t) ~ A[l - cL(t)- 2 x 2 }L e (t) F(x) (32) 

when d > d c + 2. In eqs. ( |30|) (pl|) and ( p2[) a, a' and c are constants 
and d c = 9 + 2 is a critical dimensionality playing a role similar to that 
of d c in the quench at jjl\. 

^ 3 = (r = 0,g<0,A>0) 

As in the NCOP model, the quench at /Z3 is very peculiar because the 
necessary condition @ suggests that metastable relaxation to $(f, t = 
00) = is not ruled out, but the sufficient condition (|10|) does not apply 
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here since it reduces to S(t) < which is never true for finite times. We 
will show that both the dynamics leading to stable and to metastable 
equilibrium are consistent with the model equations. 

When stable equilibrium is reached, since in this case lim^oo S(t) = 
— | 7^ from eqs. (§), (0) and (§) we observe that only lim t _ i . 0O (3(t) = 
— oo is possible. Therefore, from eq. (H), letting S(t) ~ — |, it is found: 

Q(*)~-(^*M* (33) 

and eventually: 

C(k,t) ~ Ae [{d - e)lnt] ^ x2 L e (t) F(x) (34) 

As in the $ 4 theory standard scaling is broken in the limit of large N. 



Instead, a multiscaling symmetry shows up in eq. (|34|). By considering 



the peak of the structure factor k m {t) it has been found: 

[k m {t)L{t)f ~ ~^-lnt (35) 

In other words in the late stages of the quench, two distinct lengths /c" 1 
and L(t) exist and differ only by a logarithmic factor. The standard 
scaling symmetry is broken by this feature and, instead, multiscaling 
holds @. 

If metastable equilibrium is approached, however, assuming limj^oo (3(t) 
0, from eq. (f|) and eq.(9) of I, to leading order, we find: 

Q(t) * 9 -pL(tr d ^ + c (36) 
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where a = r(g ^ +4) and c are constants. 
Computing the ratio 

f3(t) ~ L(t) e - d+2 + cL(t)- 2 (37) 

we deduce that this solution is consistent only for d > d c = 9 + 2. In this 
case we find: 

C(k, t) ~ Ae [9AaLit)dc ~ d+cL{t) ~ 2]x2 L e (t) F{x) (38) 

Hence: 

C(k, t) ~ A[l - gAaL{t) d ^ d x 2 ]L e {t) F(x) (39) 

for d c < d < d c + 2, and: 

C(jfe, t) ~ A[l - bL(t)- 2 x 2 }L e (t) F(x) (40) 

for d < d c + 2 (with 6 = c for d < d c + 2 and b = a + c for d = d c + 2). 

Metastable equilibrium is approached with the same asymptotic dy- 
namics as in a quench at //2, for d > d c . 

On the other hand, for d < d c , none of the asymptotic forms (0) and 
(H) are consistent. In fact let us try firstly with (0): 

s(t) ~ g^" (4i) 

If metastability is approached: 

Qit) ~ -bigV (42) 

in the late stage. Since: 

lim Qit) = +oo (43) 
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from eq. ( f4"2|) it is: 

lim Q(t) = -oo (44) 

Statements (^) and (|44[) cannot be true simultaneously. 

On the other hand, even the asymptotic form @ can never be 
consistent with the requirement S(t) — ► 0. In this case, in fact, eq.( 9) 
of I would read: 

Q(t) = -\g\S(t) ~ -|$||Q(t)| ^e^ 1 (45) 

asymptotically This equation, however, quickly leads to a diverging 
Q(t), which is not consistent with metastability (i.e. S(t) — > 0). 

To summarize this /13 case we observe that the situation is not qual- 
itatively different from the NCOP case, in that for d < d c stable equilib- 
rium is always achieved, while, for d > d c , both stable and metastable 
solutions are possible and one passes from the former to the latter by 
changing A, the variance of the initial condition. 

/x 4 = (r < 0,g,A>0) 

In this case statement (|D|) prevents metastability. Therefore only 
lim t ^ 00 /3(t) = —00 is consistent and, from eq. (§), letting: 

S(t) K - g + ^^ (46) 

for long times, the same form (^) is found for Q(t), as in the previous 
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case. Hence: 

C(k,t) ~ Ae^ d -^ lnt ^ x2 L e (t) F(x) 

as for stable relaxation at ^3. 



(47) 



/i5 = (r>0,g>-Vf,A>0) 



Now linii^oo S(t) = 0. Hence, solving eq. (9) of I for long times we 
find: 

Q(t) ~rt + c (48) 

where c is a constant. Therefore C(t) ~ y/rL 2 (t) is the dominant length 
and, by means of eq. fl26|) we obtain: 



C(k,t) = Ae ^)F(i) 



A 



x 



rC 2 (t) 



F(x) 



(49) 



This result shows that, differently from the corresponding NCOP case, 
scaling holds controlled by the length C(t), which grows as £«, with z = 2 
in any dimension. 

fit = (r > 0, g < -±M , A > 0) 



According to statements (0) and QXD|) if the dynamics leads asymp- 
totically the system to a state with S(t) < *A -, metastable relax- 
ation in $ = occurs. Practically this is always achievable by choosing 
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an initial condition with 5(0) sufficiently small. In this case the equa- 
tions can be solved asymptotically as in a quench at ^5. The same results 
are obtained: 



C(k,t)~A 



1 



:i l 



F(x) (50) 



rC 2 {t) 

In this case, as compared to /i3, the stronger character of metastability 
eliminates the critical dimension: metastable relaxation occurs in any 
dimension when S(0) is small. Since (0) is only a sufficient condition, 
eq. (|50|) can in principle hold asymptotically even for S(t) larger than 



|g|-yV-4Ar 

2A 

When S(0) is sufficiently large, however, stable equilibrium must be 
obtained. In this case, proceeding as for /i 3 or /z 4 , the same multiscaling 
solution (|33) is obtained in the late stages: 



C{k,t) ~ A e [(d-°)^* 2 L°(t) F(x) (51) 

2 Summary 

In this paper we have extended the solution of the spherical <J> 6 model to 
the conserved case (model B), with zero temperature. 

To summarize the results of the present paper we observe that the 
<3> 6 model shows either scaling or multiscaling in the asymptotic regime 
depending on the parameters \x of the Hamiltonian. The asymptotic 
multiscaling, peculiar to the large-iV limit, shows up, as expected, in 
those cases when stable equilibrium is reached in a potential with two 
degenerate minima (see I, fig.l). The standard scaling behaviour, on 
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the other hand, is controlled either by L(t) or by C(t), with different 
exponents and scaling functions in the two cases. 

The scaling behaviour controlled by L(t) is always induced by the 
presence of the trivial fixed point of simple diffusion at fx , which can be 
attractive on the whole r = axis. For g > this line represents the edge 
of the sector /i 5 where C(t) prevails asymptotically. As a consequence the 
trivial fixed point is attractive only above a dynamical critical dimension, 
while below the structure factor scales with C(t) as at ^5. 

The negative part of the g axis, on the contrary, is the intersection set 
of two sector (// 4 and /x 6 ) where multiscaling holds. Therefore the trivial 
fixed point (which represent now metastable relaxation) competes, now, 
with the multiscaling fixed point and, again, it can be attractive only for 
dimensions which exceed the critical one, in analogy with static critical 
phenomena. 

Finally, for what concerns the /xg region, it behaves as a transition 
sector between the /X4 and /X5 sectors of the phase diagram, in that both 
multiscaling and scaling induced by C(t) are found (the last leads to 
metastable equilibrium). Since r > the trivial fixed point is never 
attractive and no critical dimension is found. 

As regards metastability, it occurs in strict analogy with the non- 
conserved case, apart from a greater technical complexity. In both cases, 
in fact, a similar criterion links metastability to a local property (its 
curvature) of the potential around the metastable extremum. We con- 
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elude that the kind of dynamics (whether conserved or not) does not 
affect the occurrence of metastability which is due, on the contrary, to 
the analytical properties of the Hamiltonian.This is an interesting fea- 
ture which suggest an unified description of physical processes involving 
thermodynamic metastability. 
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4 Figure caption 

In figure 1 various shapes of the potential V(|0|) are schematically shown 
as a function of the parameters \x = [r, g, A], for A > 0. 
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